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I. INTRODUCTION 

Topological order is responsible for interesting new 
states of matter that do not fit into the standard 
symmetry-breaking pictured For decades, the Landau 
paradigm successfully described systems by looking at 
the underlying symmetries, with phase transitions oc- 
curring between phases with different symmetries. How- 
ever the integer quantum Hall (IQH) effect showed this 
approach to be incomplete, since it exhibits transitions 
between phases of the same symmetry. These phases are 
instead distinguished by topological order, with gaplcss 
modes localized at domain walls between regions of dif- 
ferent topological order. 

For systems with a non-zero bulk gap at all points in 
the Brillouin zone, it is possible to define a topological 
invariant of the Hamiltonian. Systems with non-trivial 
topological invariants are termed topological insulators 
(TI) and topological superconductors^— The invariants 
are robust to smooth deformations that do not close the 
bulk gap and underlie the precise quantization of re- 
sponse functions in topologically ordered systems. This 
was first realized in IQH states, where the Hall conductiv- 
ity can be expressed as the first Chern number of the (7(1) 
vector bundle of Bloch states^ - — For time-reversal invari- 
ant (TRI) systems in two-dimensions (2D), a Z2 topolog- 
ical invariant distinguishes between the vacuum (trivial 
phase) and the quantum spin Hall (QSH) stated— In 
three-dimensions (3D), there are four Z2 invariants de- 
scribing TRI systems, one of which distinguishes between 
the vacuum and a strong topological insulator and is ro- 
bust to disorder.— ~— The quantized magnetoelectric re- 
sponse may be written in terms of this Z2 invariant i 18 ' 19 

Although a topological invariant is an abstract quan- 
tity defined for a fully periodic system, it is manifested 
physically as mid-gap surface states^ In IQH systems, 
the quantized Hall conductance can be formulated in 
terms of the the number of chiral edge states! 21 ' 22 Sim- 
ilarly, the 2D/3D Z2 invariant for TRI systems deter- 
mines whether there are an odd or even number of helical 
modes/Dirac cones at a given edge or surfaced In the 



cases above, the edge states smoothly connect the bulk 
valence and conduction bands and the number of such 
modes is protected by the topological invariant: they can- 
not be deformed into a single bulk band unless the bulk 
gap closes. By contrast, edge modes in an ordinary sys- 
tem do not traverse the bulk gap and are thus susceptible 
to localization by disorder. For superconducting systems, 
the topological invariants determine the number of Majo- 
rana modes localized at the edge or in vortices i 24 i 25 These 
states are at zero energy and are protected by particle- 
hole symmetry and the superconducting gap. Systems 
such as p + ip superconductors (SC) in the "weak pairing 
phase" or SC-TI hetero junctions can support Majorana 
modes which obey non-Abelian statistics i2£r— 

The goal of this paper is to derive a rigorous connection 
between the bulk invariants and the surface dispersion. 
A heuristic way to understand this bulk-boundary corre- 
spondence is as follows. Consider a domain wall between 
two bulk insulators with suitably defined topological in- 
variants that take the values vl and vr ^= vl in some 
regions. Since the value of the invariant cannot change 
for finite energy gap, this means the bulk gap closes at 
some interface. Mid-gap excitations can thus exist, but 
they are confined to the interface by the bulk gap in the 
other regions. This argument applies to domain walls 
between regions with different values of the invariant, of 
which an edge is a special case where one of the regions 
is the vacuum (trivial phase) i 29 i 30 

In light of recent interest in topological insulators and 
superconductors, it would be useful to formalize the re- 
lation between bulk topological quantities and proper- 
ties of mid-gap edge states. This connection has been 
proven specifically for IQH states on a square lattice by 
deriving a winding number for the edge states^ Another 
approach using twisted boundary conditions has the ad- 
vantage of including interactions and disorder, but can- 
not prove that the states exist at an open boundary^ 

There has also been recent progress on analytic so- 
lutions of edge states in topological insulators^— and 
topological superconductors.— These calculations are of- 
ten based on models using a specific Dirac Hamilto- 
nian. Dirac systems are ubiquitous in condensed mat- 
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ter and particle physics systems and give rise to many 
exotic states. For example, every single-particle topo- 
logically ordered system can be realized with a Dirac 
Hamiltonian4»22 They are used to model a variety of 
systems including polyacetylene, quantum Hall insula- 
tors, graphene, topological insulators and superconduc- 
tors, etctfWHttZ 

In this paper we deepen the understanding of surface 
states by deriving their dispersion, effective theory, and 
chiral properties. Our work applies specifically to tight- 
binding Dirac Hamiltonians with nearest-layer interac- 
tion. For these systems we present a prescription for the 
edge states spectrum and prove the bulk-boundary cor- 
respondence. In addition, we derive a simple geometric 
method to calculate the energies and penetration depth 
of the edge states analytically. 

The organization of the paper is as follows. In sec- 
tion UU we introduce the bulk quantities of a lattice 
Hamiltonian that determine topological behavior. In sec- 
tion IIII1 we state and prove the two main results of the 
paper, Theorem 1 relating the parameters of the bulk 
Hamiltonian to the surface spectrum in a geometric way, 
and Theorem 2 proving the bulk-boundary correspon- 
dence between chiral edge states and the Chern number. 
In section [TV] we demonstrate the range of applicabil- 
ity of our theorems and give examples of topologically 
ordered systems. We also show how the bulk Z2 invari- 
ant for a time-reversal symmetric insulator relates to the 
number of surface Dirac cones in example IIV CI In sec- 
tion [Vj we extend the results from lattice Hamiltonians 
to continuum quadratic Hamiltonians, with discussions 
on its implications. In closing, we discuss the possible 
extension of the work to other classes of topological su- 
perconductors beyond IQH and TRI systems. 



II. CHARACTERIZATION OF THE 
NEAREST-LAYER HAMILTONIAN 

To study a system with edges, consider a 2D/3D crys- 
tal that terminates on a line/surface. Translational sym- 
metry is thus broken in the direction normal to the edges. 
However, we assume it is unbroken parallel to the surface, 
and the corresponding momentum k|| is a good quantum 
number. In this way, any higher dimensional system can 
be decoupled into a family of one-dimension (ID) prob- 
lems parameterized by kii . 

The Dirac Hamiltonian in momentum space H (k) can 
always be expressed as a linear combination of gamma 
matrices, H (k) = h(k) • T. Here T is a vector of the her- 
mitian gamma matrices (independent of k) which satisfy 
the Clifford algebra TT^ + T^T 1 = 25 ij . h is a real vec- 
tor that maps the Brillouin zone to a closed curve in 
a g-component vector space, where there are g gamma 
matrices T % . The Pauli matrices are examples of gamma 
matrices: any 2x2 traceless matrix can written as a h ■ cr 
where h is a 3-component vector and <x = (a x , a y ,cr z ). 

Squaring the Hamiltonian gives H 2 = (h ■ T) 2 = |h| 2 . 



The eigenvalues of H, given by E(k) = ±|h(k)|, can be 
thought of as the distance of the vector h to the origin. 
If H(k) describes a band insulator with a bulk gap, then 
the locus of points traced by h(k) never intersects the 
origin. 

Let us consider Dirac Hamiltonians with coupling be- 
tween neighboring layers: 

n = E <k„r- [b*n-i,k„ +b°* n>k|I +b** n+llk „] , 

n,k|| 

(1) 

where n labels the layers. Both b and b° are depen- 
dent on k|| but we will not write this dependence explic- 
itly. 'I'n.kii is a vector of quasiparticle annihilation oper- 
ators at layer n that captures all the degrees of freedom 
(i.e. spin, pscudospin) at every site. Fourier transform- 
ing (n — > fcj_) in the direction away from the edge, the 
bulk Hamiltonian becomes 

U = Y, *k [h(k) ■ r] tt k , (2) 

fc_L,k|| 

with 

h(k) = be-' lk± + b° + b*e ik± 

= b + 2b r cos/c ± + 2b l sin/c ± , (3) 

where b r and b 1 are the real and imaginary components 
of the vector b respectively. We point out that b and b° 
are independent of fcj_ 

The curve traced out by h(k) for fixed kii is an ellipse 
living in the plane spanned by b r and b\ b° can be 
decomposed into a component b*]_ that is normal to this 
plane, and bjj that lies within it. bjj shifts the ellipse 
within the plane, while b 1 ^ lifts it uniformly. It will be 
useful to define 

h,| (k) = bjj + 2b r cos k± + 2W sin k± (4) 

as the projection of h(k) on to the 2D plane spanned 
by b r and b\ Note that this plane contains the origin, 
while the plane containing h is offset from the origin by 
the vector bj_. Since h|| maps the Brillouin zone to closed 
curves, it can be divided into two classes: ellipses that 
enclose the origin, and ellipses that do not. 

III. EDGE STATE AND BULK-BOUNDARY 
THEOREMS 

As we shall prove in this section, the behavior of h(kj_) 
completely determines the topological nature of the sys- 
tem and holds the key to understanding the relation be- 
tween existence of edge states and bulk topological in- 
variants. This section contains the main result of the 
paper, where we prove two theorems, one relating the 
spectrum of edge states to h, the other connecting h to 
a bulk topological invariant. 
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FIG. 1. An illustration of Theorem 1. The gray ellipse is 
traced out by h(fcj_) = b° + 2b r cos k± + 2b ! sin k± for a fixed 
parallel momentum ky [Eq. |(3J]. The dotted ellipse (hy) is h 
projected on to the plane spanned by b r and b 1 . The displace- 
ment of the ellipse h from the dotted ellipse hii is given by bx, 
the component of b° perpendicular to this plane. Theorem 1 
says that an edge state exists if and only if the dotted ellipse 
encloses the origin (which holds true for the diagram above), 
and its energy is determined by the displacement \h° ± \. 

A. Edge state energy 

Theorem la. The system has mid-gap edge states if 
and only ifh^(k±) encloses the origin. 

Theorem lb. The energies of the edge states are 
given by the distance from the origin to the plane contain- 
ing h, i.e. E s = ±|bj_|. When the gamma matrices are 
the Pauli matrices, the energy of the left edge state (semi- 
infinite slab with n > 0) is given by: E s = b° ■ i^f^ti • 

Here we only sketch the main ideas behind two equiva- 
lent proofs of the theorem, leaving the full details to the 
Appendices. We present two approaches to this problem: 
a proof utilizing Green's functions^ (c.f. App. [SJ and 
a proof via transfer matrices^ {c.f. App. [B]). In this 
section, we consider one block of the decoupled system 
corresponding to fixed k|| . 

We begin by writing the Green's function of the sys- 
tem without edges, where the full translational invari- 
ance makes a momentum space representation possible. 
A system with edges is then created from the fully peri- 
odic system by deleting the couplings between one pair 
of neighboring sites. The poles in the Green's function 
G(E) at mid-gap energies E indicate the presence of edge 
states. 

The bulk Green's function is given by 

Co(^)=£Mf> (5) 

i 

where i sums over the energy eigenstates of H(kj_). Since 
we are interested in a boundary localized in real space, it 
is necessary to Fourier transform the bulk Green's func- 
tion. For a system of size L this results in a L x L 
block matrix Go(E; yi, yj), where each block corresponds 
to mixing between the lattice sites yi and yj. 

Next we write an expression V for the terms in the 
Hamiltonian that create the boundary by subtracting the 
hopping terms between sites y\ e f t and bright- For mod- 
els with nearest neighbor interactions, the only non-zero 
matrix elements of V are those between yi e f t and ?Mght- 



The Dyson equation gives an exact expression for the 
open boundary Green's function G in terms of the bulk 
Green's function Go and the cuts V needed to take the 
system from one geometry to the other: 

G(E) = (l-G (E)V)~ 1 Go(E). (6) 

The poles of G(E) occur when the edge state energy sat- 
isfies Dct [I — Gq(E)V] — 0. If an edge state exists with 
wavefunction u, it must satisfy (/ — GqV)u = 0. We 
note that this is the same as the Schrodinger equation 
(E — Hq — V)u = 0. The benefit of the Green's function 
formalism is that it reduces the problem to only edge de- 
grees of freedom and enables an analytic solution. This 
implies the following two statements: 

J dk± ^^°' (7a) 
r , h, • b* 

J dk±e ih^r^ 71 - (7b) 

These conditions are satisfied if and only if hy encloses 
the origin, and the edge mode energy is given by ±|bj_|, 
where the sign is given by the orientation of hy . 

To prove Thm. 1 using transfer matrices, we con- 
sider a semi-infinite system with unit cells labeled by 
n = 1,2,3, etc. We seek a solution ip n to the single- 
particle Schrodinger equation: 

b • r^n-i + b° • Y^ n + b* • I>„+i = Ei'n (8) 

for n > 1. At the edge site n = 1, we have b° ■ Tipi + 
b* • r-02 = Eipi. This condition is enforced by applying 
Eq. JSJ for n = 1 but stipulating that ip = 0. The re- 
cursion relation ([8]) relates ip n +i to ip n and tp n -i. Hence 
given ipi (and ipQ = 0), we can recursively calculate all of 
ip n and construct the wavefunction. 

An edge state requires ip n to be exponentially decaying 
as n increases, hence the solution ip takes the form: 

4> n =u a X2+u b X^ (9) 

where u a = — u bl and A Q , are complex with |A a |, |A(,| < 
1. Algebraically, this is equivalent to having E = ±|b^| 
and finding two roots within the unit circle of the func- 
tions L(X) or L(A), defined as 

L(A)=h )1 HlnA)-(v 1 +zv 2 ), 
L{\) = h||(-ilnA) • (vi -iva), 

where vi , V2 are two orthonormal vectors that live in the 
plane of h||. When A = e lk± lies on the unit circle, L(X) 
and L(X) trace out the ellipse h(fcjj in the complex plane 
clockwise and counterclockwise, respectively. Because of 
this property, the number of times h(fcx) wraps the origin 
determines the number of zeroes of L(X) and whether the 
two solutions X a ,b in © exists. In Appendix [Bl we pro- 
vide the full details bridging these steps, and also com- 
pute the sign of the edge state energy as well as their 
penetration depth. 
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B. Bulk Chern number and chiral edge 
correspondence 

In this section we prove Theorem 2, relating the bulk 
Chern number v with the number of chiral edge modes 
for 2 x 2 Hamiltonians. 

Theorem 2. A chiral edge mode exists for a 2D bulk 
insulator if the bulk has a non-zero Chern number, i.e. 
h(k) wraps the origin. The number of chiral edge modes, 
counterclockwise minus clockwise, is given precisely by 
the Chern number. 

When the irreducible representation of T are 4 x 4 or 
larger, it can be shown that the Chern number is always 
zero. The edge states of any surface always appear in 
pairs with energy +E S and — E s and so the number of 
clockwise and counterclockwise chiral modes are always 
equal. We are particularly interested in 2 x 2 Hamiltoni- 
ans because they can have nonzero Chern numbers and 
support chiral modes. 

Consider an insulator in two dimensions whose Hamil- 
tonian is written as a 2 x 2 traceless matrix: H(k x , k y ) = 
h(k x ,k y ) ■ cr. Because the bulk gap of an insulator is 
non-zero, h is non-zero at all points in the Brillouin zone. 
Hence H (k) is a map from the Brillouin zone (torus) to a 
set of non-zero vectors with 3 components (K 3 — {0}), and 
such maps can be characterized by a v £ Z topological 
invariant, known as the Chern number 4£ Hamiltonians 
with different Chern numbers v cannot be deformed into 
one another without closing the bulk gap. In this con- 
text, the invariant v determines the number of times the 
torus h(k) wraps around the origin. 

To examine the edge states for an arbitrary edge, say 
one parallel to y, we analyze the spectrum as a function 
of feii = k y . The torus h(k) can be divided into a family 
of loops h(k x )\k , each at a fixed value of k y and giving 
information of the edge state at that momentum. 

Before proceeding to the technical proof, we present a 
geometric argument with the aid of Fig. [5J which shows 
an example of a bulk insulator with Chern number v = 1. 
The important loops of fixed k y are highlighted in black. 
Since v is nonzero and the torus wraps the origin, it is 
always possible to find two loops that are coplanar with 
the origin, one of which encloses the origin and one that 
does not. In this example, the latter case occurs at k y = 
0, indicating no mid-gap edge states at this k y . As we 
scan through different values of k y , the loop moves out 
of this plane. At some critical momentum k c (given by 
■| in Fig. [2J , the projection of the loop onto this plane 
intersects the origin and an edge state emerges from the 
bulk conduction bands. At k y = tt, the loop is coplanar 
with the origin and encloses the origin, indicating zero- 
energy edge states at this value of k y . As the plane of 
the loop passes through the origin, the energy of the edge 
state changes sign. The presence of edge modes for this 
range of momentum is shown as orange shading in Fig. [2] 
Eventually at some critical momentum k v (given here 
by ^), the loop moves away from the origin and the 
edge state disappears in to the bulk valence band. Since 




(a) 




(b) 



FIG. 2. (Color online) Illustration of bulk-boundary cor- 
respondence. Figure 2(a) shows the torus traced out by 
h(k x ,k y ) for a bulk insulator with Chern number v = 1. 
Each black loop maps out h(k x )\k y for fixed values of k y , 
the thick black lines guiding the eye to important loops. Set- 
ting k y = n gives the black loop on the right that encloses 
the origin, meaning there are zero-energy edge modes at this 
value of k y . At k y = 0, the black loop on the left lies in the 
plane of the origin without containing it, indicating no edge 
mode at k v — 0. The black loops on the top and bottom 
|) have projections which intersect the origin, 



{kv — g , k, 

indicating the values of k y where the edge band merges with 
the bulk bands. Figure [2(b)| shows the band structure of the 
system with the edge mode drawn in orange. The model pre- 
sented here is ap+ip superconductor described in section ll V Bl 
[see Eq. (|17[) ] with parameters: t = 1, Ao = 3, y, = 1. 



the edge state energies at k v and k c have opposite signs, 
the edge band connects the bulk valence and the bulk 
conduction bands. 

Formally, we can describe each loop h(k x )\k y by the 
Berry phase <f>(k y ) living in a circle [0, 2ir] with ~ 27r4£ 
The Berry phase can be formulated in various ways: 



p2tT 

<t>{ky) — - j> dk x A x (k x ,k y ) 



dk'i dk x F(k x ,k') 



(lla) 

(lib) 
(11c) 



where Aj(k) = i(ipk\djipk) is the Berry connection of the 
filled energy states of H(k), F = d x A y — d y A x is the 
Berry curvature. Geometrically, <f> is half the oriented 
solid angle f2(h) subtended by the loop h(k x ) as seen 
from the origin. The integral of over the entire Bril- 
louin zone gives the Chern number: § m F = v. Both 
4> and k y live on a circle, and from Eq. (jllbp . 0(fc y ) is a 
map S 1 y S 1 with winding number v. 

At the values of <fr(k y ) = or ir, the origin is in the 
plane of the ellipse h(k x )\k y , and lies outside or inside the 
ellipse respectively. Hence there is a zero energy edge 



state when 4>(k y 



and no edge state if <f>(k y ) = 



(or 27r). The family of loops as k y is varied connects 
these two special cases smoothly. For example, the upper 
critical momentum k c has Berry phase < <p{k c ) < it, 
while the lower critical momentum k v has Berry phase 
it < 4>(k v ) < 2ir. Thm. lb says that if an edge state 
exists, < < 7r implies it has energy E s > 0, and 
it < <f> < 2tt implies E s < 0. Therefore in between k c < 
k y < k v , a gapless (counterclockwise) chiral mode must 
exist connecting the bulk bands. 

For an insulator with Chern number v, the Berry phase 
4>(k y ) goes from to 2ttv as k y is varied from to 
27r. Each time the phase <j>{k y ) winds around the circle, 
27ra — > 27r(a + 1), a pair of critical momenta (k ca ,k va ) 
defines a range in which a chiral mode connects the bulk 
valence and conduction band, k ca < k y < k va . This 
chiral mode is counterclockwise as the phase <p increases 
by 27T. Similarly, there is a clockwise chiral mode as <j> 
decreases by 27r. Therefore, the total number of chiral 
edge modes (counterclockwise — clockwise) is given by 
the Chern number of the bulk Hamiltonian. 



C. Discussion 

Theorem 1 gives a simple way to compute the spectrum 
of edge states from properties of the bulk Hamiltonian. 
The existence of zero-energy edge states is determined 
by whether or not the ellipse traced by h|| encloses the 
origin. Intuitively, the size of the ellipse is a measure 
of the coupling strength b between neighboring layers, 
while the in-plane displacement of the ellipse bjj is a 
measure of coupling within the layers. From this per- 
spective, Thm. la says that an edge state exists if the 
nearest-layer coupling is 'stronger' than the intra-layer 
coupling. This is a straightforward extension of the edge 
states of polyacetylene, a ID chain with alternating bond 
strengths t ^ t' , which supports an edge state if the chain 
terminates on the weaker bondj^ 

The argument presented above can also be used to cal- 
culate the penetration depth £ of the surface states: 



21n(l/|A|) 



(12) 



a is the distance between layers, A is the characteristic 
decay parameter such that the wavefunction decays as 
ipn ~ A™ in the bulk. |A| is the larger of |A a |, |A{,| [defined 




FIG. 3. Determining the penetration depth from the ellipse 
h||. The distances from the foci of the ellipse h||(fcx) to the 
origin determine the characteristic decay parameter A, which 
in turns gives the penetration depth f = — o/21n|A|. In the 
case where hy traces a circle, |A| = d/r, where d is the dis- 
tance of the origin to the center of the circle, and r is the 
radius of the circle. In the general case hy traces an ellipse, 



|A| 



, where M and m are the major and minor 
diameters, / is the distance between the foci IF1F2I, and I is 
the sum of distances \OF\\ + |OF 2 |. 



in Eq. ([9])]. |A| is always less than one and is determined 
by the location of the origin inside the ellipse h||(fcj_). 
When the origin touches the edge of the ellipse, A has 
unit modulus and £ tends to infinity, indicating a bulk 
propagating mode. At this point the surface spectrum 
ends and merges with the bulk bands. The decay param- 
eter can be computed by 



|A| 



M + m 



(13) 



where M and m are the major and minor diameters of 
h||(fcj_) respectively, / = V M 2 — m 2 is the distance be- 
tween the foci of the ellipse, and I is the sum of the dis- 
tances from the origin to the two foci. (See Fig. [3J) 

The bulk-boundary correspondence described in 
Thm. 2 holds even in the presence of weak interactions. 
The number of edge states is a topological property of 
the bulk and is robust to small perturbations that do 
not close the bulk insulating gap. Such weak interac- 
tions have the effect of altering the edge state dispersion 
from the simple formula derived in Thm. 1. Neverthe- 
less, we may still view the result of Thm. 1 as a starting 
point where the weak interactions are added perturba- 
tively. With strong interactions, the electronic excita- 
tions may no longer behave as a Fermi liquid, leading 
to a break down of the bulk-boundary correspondence 
described here. For example, electron-electron interac- 
tions may drive the ID edge states of QSH insulators 
to become Luttinger liquids, which have been shown to 
break the topological distinction between the insulating 
and helical edge states i 47 ' 48 



IV. APPLICATIONS OF THEOREM 1 

In this section, we illustrate how Theorem 1 may be 
used the edge state dispersion in various systems. 
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A. Example: graphene 

As an illustration of how Thm. 1 can be used to cal- 
culate edge state energies, we examine the zigzag edge of 
graphene. Because the bulk energy bands are degenerate 
at two inequivalent points in the Brillouin zone, there 
is no well-defined topological invariant and Thm. 2 does 
not apply. 

Consider the tight-binding model for graphene on a 
honeycomb lattice with primitive translational vectors 
ai and a2 taken to be 120° apart. The zigzag bound- 
ary parallel to either ai or a2 is known to support edge 
modes, while the armchair boundary parallel to ai — &2 
has no edge modes. Using the reciprocal space coordi- 
nates fci = k • ai , &2 = k • a2 , the tight-binding Hamilto- 
nian takes the form 



H = t 



1 





-iki 



i 



-iko 



ik-2 







(14) 



written in the (A, B) basis, where A and B are the in- 
equivalent sublattices. The Hamiltonian can be written 
in terms of the Pauli matrices as H = h • cr with 

h(k) = i(l + cos fci + cos It2, — sin fci + sin fc'2, 0). (15) 

For a zigzag edge parallel to ai, we examine the curves 
traced by h for fixed fci 4^ The /^-dependent terms in h 
are written as 

2b r cos fc2 = (t, 0, 0) cosfc2, 
2b l sinfc 2 = (0,t,0)sinfc 2 , 
while the remaining fc2-indepcndcnt terms become 
1 '-fcosfci,-tsinfci,0), b°=0. 



bfl = (t 



h|| lies in the plane spanned by b r and b\ and for this 
system is a circle in the xy-plane. It encloses the origin 
only if the magnitude of b° is less than t, i.e. when 



|1 + e lkl | < 1 or equivalently when 2(1 + cos fci) < 1. 
Thus for 4p < fci < 4^, h|| encloses the origin and the 
system has mid-gap edge modes. According to Thm. lb, 
the energies of these edge modes are given by ±|bj_|, 
which is zero. We have shown that the zigzag boundary 
supports zero energy edge states when ^ < fc x < 4^, in 
agreement with previous work i 41 ' 50 

To see why such edge modes do not exist at armchair 
boundaries, consider the edge parallel to ai — a2. By 
rewriting the Hamiltonian in terms of fc|| = fci — fc2 and 
kx = fci, we can again examine the behavior of hn at 
fixed fc|| . In this case, it is possible to show that h|| never 
encloses the origin, and thus never satisfies the condition 
for mid-gap edge states. 



B. Example: p + ip superconductor 

We study a p + ip system, which in the weak-pairing 
phase is characterized by chiral Majorana modes at the 



edgei^i We consider a model on a square lattice with p- 
wave pairing between nonrelativistic spinless electrons, 
shown in Fig. [5] The bulk Hamiltonian is given by 



H(k) = 



-6c- 



A, 



(16) 



where H is written in the (ck,cl k ) T basis. The kinetic 
energy is given by £k = — 2i(cos k x + cos k y ) and the pair 
potential is given by Ak = Ao(sin k y + i s'mk x ), with Ao 
real. Eq. (fl6|) can be expressed as H = h • r, where t % 
are the Pauli matrices in Bogoliubov-de Gennes (BdG) 
space: 

= (Ao sinfcy, — Ao sinfcj;, — 2t cosk x — 2tcosk y — /i). 

(17) 

For a system with edges parallel to x, we examine the 
ellipses traced by h(fcj,) for fixed values of k x . Following 
Eq. ([3]) , the fc y -dependence of h is decomposed as 

2b r cosfc y = (0,0, -2t) cosk y , 
2b' sinfc y = (A , 0, 0) sinfcy. 
The remaining fc. y -indcpendent term is decomposed as 
bjj = (0,0,-2tcosk x -n), 
b° ± = (C-Aosinfc^O). 

For this model h||(fc y ), which lies the plane spanned by 
b r and h l , is an ellipse in the ccz-plane. The condition 
for the ellipse to enclose the origin is \2t cos k x < \2t\ , 
which is only possible for the range of chemical potential 
< 4\t\. According to Theorem 1, when this condition 
is satisfied, the system has mid-gap states at the left edge 
with energy 



E s = b° ± 



b r x b i 



b ' r x b * 
A sin k x , 



K ■ (-v) 



assuming tAo > 0. The right edge state energy is given 
by — E s = —Ao sinfcy. The two edge states can become 
degenerate at E s = at either k x = or 7r for an appro- 
priate range of [i: the degeneracy occurs at k x = when 
-4 < fi/t < 0, and at k x = n when < fx/t < 4. The 
bulk gap closes and there is a transition at fx/t = or 
±4. 



C. Example: 3D topological insulator 

In this section we show that the surface states of a 
strong topological insulator have an odd number of Dirac 
cones, and derive an effective surface theory. As an ex- 
ample of a topological insulator, we use a model on a 
cubic latticed 



H 



VT 



sin fc, 



M 



cos fc^r* (18) 
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for fx = x, y, z. In the basis of Dirac matrices T = 
(t z u x ,r z a y ,r z a z ,t x ), we can write H = h • T, with 

h(k) = (v s'm k xi v sin k y ,v sin k z ,m()s.)), (19) 

where the "mass" is given by m(k) = M — i^cos/c^. 
(The 'fifth' Dirac matrix is t v , but it is never used in 
this model.) For the (001) surface parallel to x and y, we 
solve for the edge states as a function of k x ,k y : 

h(fc z )| fe k = (v sin k x , v sin k y , 0, M — t cos k x — t cos k v ) 

+ (0,0,0, -t)cosk z + (0,0,i;,0)sinfc a , 

which lies on the 34 -plane (in T space) displaced by 
w(sin k x , sin ky) from the origin. An edge state exists if 
and only if \M — tcosk x — tcosfcjj < |t| with Dirac cone 

spectrum E s = ±v^J sin 2 k x + sin 2 k y . 

At any of the surface time-reversal invariant momenta 
(TRIM) (k x ,k y ), the ellipse h|| (fc z )|fc xi fc H has one of its 
diameters stretching from m(k x ,k y ,0) to m(k x , k y , tt) 
along the T 4 axis. This ellipse encloses the origin 
if and only if the two endpoints straddles the ori- 
gin; equivalently, a Dirac cone appears at the surface 
TRIM if the mass in the bulk Hamiltonian changes sign: 
m(k x , ky, 0)m(k x ,k y ,ir) < 

In this basis the time-reversal operator has the form 
= —ia v K, where IC is the complex conjugation opera- 
tor. The spatial inversion operator is II = T 4 = t x . At 
the eight bulk TRIM, the Hamiltonian commutes with 
the inversion operator: [H ( k^RiM ) > II] = 0. In fact, H is 
a multiple of II: 



-ff(k T RiM) = m(k T RiM)n. 



(20) 



Hence, the sign of the mass determines the parity eigen- 
value of the pair of filled states. The strong Z2 topo- 
logical invariant is the product of the parity eigenvalues 
at these TRIM, and hence v$ = Htrim sgnm(k)j^ It is 
clear that the bulk Z2 invariant dictates whether there 
are an even or odd number of Dirac cones on the sur- 
face^ 2 . 

Generically, all TRI topological insulators written in 
4x4 Dirac matrices take on a similar form to (fT5|) . 
and most importantly they satisfy (f2T)|) at the eight bulk 
TRIM. Hence it is straightforward to establish the bulk- 
boundary correspondence for 3D TRI topological insula- 
tors. 

We now proceed to examine the surface eigenstatcs 
following the derivation from Sec. IB 71 When 1 < m/t < 
3, the model is in the strong topological insulating phase 
with a single Dirac cone at (k x , k y ) = (0,0). The surface 
Hamiltonian is of the form H s = (b^ • T)P S where P 8 
is the projector of the two surface states: P s = ^(1 — 
i(— T 4 )r 3 ) = i(l + r v a z ). To examine the Hamiltonian 
in the reduced vector space, it is useful to apply a unitary 
transformation H = UHU^ 1 where U = cxp(— ijT x ). 
The projector becomes diagonal in the new basis 



Ps 



1 



projecting onto the first and last row. The "edge energy" 
term h°, ■ T transforms as 



Vl = t7(b° • T)U~ 



—vT y (a x sin k x + a v sin k y ) 
-vT y {a x k x + a y k y ). 



Since P s is diagonal with entries (1, 0, 0, 1), it suffices to 
examine only the four corners of V®. The effective surface 
Hamiltonian can be computed)^ 



h, = v?p. 







v(cr x k) • z. 



(21) 



The basis of the surface Hamiltonian H s is, in terms of 
the basis of H, (1, 0, i, Q) T /V2 and (0, 1, 0, i) T / V2 (The 
first and last columns of ) . The spin degrees of freedom 
(cr) and the orbital degree of freedom (r) are entangled 
in the surface states. 



V. CONTINUUM HAMILTONIANS 
QUADRATIC IN MOMENTUM 

The bulk-edge correspondence stated in Sec. |H] may 
also be extended to Hamiltonians in the continuum. 
Given a translationally invariant Dirac Hamiltonian 
quadratic in momentum p = — iV, of the form: 

H(p;p {l ) = C°(p||) + C 1 (p ll )p + C 2 (p l[ )p 2 

= [c°(p || ) + c 1 (p || )p + c 2 (p||)p 2 ].r (22) 

where py and p are, respectively, the momentum paral- 
lel and perpendicular to the edge/surface. For a fixed 
momentum py, the vector 



Hp)\ 



c l p 



2 2 
c p 



(23) 



traces a parabola. This parabola lies on some 2D plane 
spanned by c 1 and c 2 , and we can always decompose c° 
and h into in-plane and out-of-plane components: c° = 
c 1 ^ + cjj and h(p) = c° ± + h|| (p). In this section we state 
the main theorem: 

Theorem 3a. An edge state (with zero Dirichlet or 
Neumann boundary condition) exists if and only if the 
origin is within the concave side of the parabola h||(p). 

Theorem 3b. The energy of the edge state is given by 
the distance of the plane to the origin, i.e. E s = ±|cj_|. 
When the gamma matrices are the Pauli matrices, the 
left edge (semi-infinite slab with x > 0) energy is given 
by: E s = c° • 



c xc 
cixc^l 



UP S U- L = -(1+rV), 



The proof of Theorem 3 is given in Appendix [C] 



A. Discussion 

The proof of Thru. 3 uses the ansatz ip(x) = u a e tKaX + 
Ube %KhX and derives the condition when Im n a & > for 
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edge states. Physically, £ = (ImAt) -1 gives the pene- 
tration depth of the edge modes. We may think of the 
continuum as a limiting case of the lattice as the lattice 
spacing a goes to zero. The ellipse hy becomes a parabola 
for vanishing a as the quantities c° = b°+2b r , c 1 = 2b l a, 
c 2 = — b r a 2 are held constant. Theorem 2 also extends 
to the continuum case as the paraboloid h(k x , k y ) deter- 
mines the number of chiral edge states. 

In addition to Dirichlet and Neumann boundary con- 
ditions, there is also a mixed type with ^'(0) = ^(O) 
for a positive number rj. The origin of this boundary 
condition comes from requiring the wavefunction outside 



(x < 0) to satisfy H v 



W + f-, where W > E 3 is 



the work function. Our analysis and result holds even 
for this boundary condition. 

Notice that the quadratic term C 2 p 2 is crucial for the 
existence of edge states. Without it, the polynomial (|C2[) 
will be quadratic and there can only be one solution for 
K in the upper half plane. The form of the solution 
-0 = ue lKX makes it impossible to satisfy either type of 
boundary condition. 

In the lattice model, we can compute the edge spec- 
trum only for certain surfaces because of the nearest- 
layer requirement. For example, we can only compute 
the {100}, {110} and {111} surface dispersion of the TI 
model on a cubic lattice. In the continuum case with 
the Hamiltonian bilinear in momentum, any surface cut 
will still yield a Hamiltonian quadratic in p±. Conse- 
quently, we can compute the surface excitation spectrum 
and wavefunctions of the system for all linear surfaces. 



B. Example: p + ip superconductor 

We use the simplest model of a p+ip superconductor >2i 

ip v ) 



2m* 
Aq(Px 



P- A (p x 

iPy) (J, - 



H(Px,Py) 

Equivalently, with H(p) = h(p) • t, 



p 

2m- 



(24) 



h (p) = ( &0Px, A Py, - \l 



(25) 



This model is isotropic, and without loss of generality, 
we take a semi-infinite plane x > with y parallel to the 
edge. At a fixed p y , the h vector becomes: 



0,A o p„, ^ -ii j 
(A ,0,0) P:c + (0,0, vl 



which lies in the xz-plane with y = Aop y . The parabola 
is concave towards the +z direction, and hence an edge 
state exists if and only if p 2 /2m* — fi < 0. Edge states 
can only exist when fi is positive, or in other words in 
the 'weak pairing phase' of p + ip superconductors. 



The edge state energy dispersion is given by Thm. 3b. 

i ? 

n C 1 X & 



E s 



AoPyV ■ (x X z) 



-A 



OPy 



Referring to Eq. (|B17j) . the edge state wavefunction is 
given by the projector P s = |(1 — it x t z ) = |(1 — t v ). 
Hence the edge states parallel to the y-axis are eigen- 
states of t v . 



VI. OUTLOOK 

In this paper, we provided two main results. Theo- 
rem 1 gives a general prescription for finding edge states 
of Dirac Hamiltonians (with nearest-layer coupling) on 
a lattice. Its range of applicability includes "accidental 
edge states" which may not be topologically protected, 
such as in boron nitride. Theorem 2 relates the bulk 
Chcrn number of a 2D insulator to the number of chiral 
edge modes. This establishes the bulk-boundary corre- 
spondence for a class of quantum Hall insulators. 

For 3D time-reversal invariant insulators, we demon- 
strated in Sec. II V Cl how the Z2 strong topological invari- 
ant determines whether there is an odd or even number of 
Dirac cones in the surface spectrum. Although we have 
used a specific TI model in the example, the argument is 
easily generalizablc for all TRI Dirac Hamiltonians. 

This work can be extended beyond quantum Hall insu- 
lators (class A) and TRI topological insulators (class All) 
to other insulators within the Altand and Zirnbauer clas- 
sification^^ The periodic table of topological insulators 
and superconductors provide an exhaustive topological 
classification of non-interacting electronic systems^^^ 
As there are model Dirac Hamiltonian o 12 ' 37 i 39 ' 51 in each 
class, our work provides the machinery to relate the bulk 
topological invariants^ to the surface propertie d 3 ' 20 of 
these systems. 
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Appendix A: Proof by Green's functions 

It is difficult to study a system with edges because of 
the broken translational invariance. We begin by writing 
the Green's function for an easier problem: a periodic 
system with no boundary— The full translational invari- 
ance allows us to work in momentum space, reducing 
the dimension of the Hamiltonian in momentum space. 
Next, the geometry of the system is changed from a pe- 
riodic to an open system by subtracting all interactions 
between a particular pair of nearest neighbors. We use 
the Dyson equation to calculate the Green's function for 
the open geometry and show that there are poles - and 
thus bound states - at mid-gap energies. The form of the 
potential required to cut the periodic system greatly re- 
duces the degrees of freedom in the problem and enables 
an analytic solution. 

We prove Theorem 1 for the case of 2 x 2 Dirac Hamilto- 
nians, which can be decomposed in terms of the Pauli ma- 
trices a*. For higher dimensional Hamiltonians, it is al- 
ways possible to find three gamma matrices T l that whose 
sub-blocks are the Pauli matrices, e.g. r z a x ,r z a y ,r z a z . 
A suitable unitary transformation rotates the Hamilto- 
nian so that it is a linear combination of these three 
gamma matrices, and the arguments of this section apply 
to each sub-block. 



1. Bulk Green's function 

We prove Theorem 1 for 2x2 Hamiltonians, which can 
be expressed in the following form: 



H{k) = E k 



cos tt sin t) e 
sin 9 e 1 ^ — cos ( 



(Al) 



Decomposing H into the Pauli matrices a 1 gives H = h er 
for 



h = Ek (sin 9 cos </>, sin 9 sin <f>, cos 9) . 
The eigenstates are given by 



KM = 



sin ; 



= (sinf h 



(A2) 



(A3) 



To express h in terms of the parameters defined in Sec- 
tion HH we assume b r and b l lie along +x and +y re- 
spectively, with b^J_ along z. The Hamiltonian takes the 
form 

Ek sin 9 cos (f> x = 2b r cos k + (bjj ■ x) x, 
E k sin 9 sin <f> y = 2b* sin k + (bjj • y) y, (A4) 
E k cos9z = b° ± . 

Here h|| lies in the xy-plane. For the remainder of this 
section, we will use the notation b r = \h r \ and hi = |b*|. 



Note that any Hamiltonian may be brought to this form 
by an appropriate gauge transformation. 

The first step in calculating the edge mode energies is 
to write the Green's function for the fully periodic sys- 
tem. Utilizing the full translational invariance, we work 
for now in the momentum representation. The a-j3 ma- 
trix element of the bulk Green's function is given by 



Gf (E;k)=J2 



E - E; 



1/3 



(A5) 



where i sums over the energy eigenstates of H^. The four 
matrix elements are given by 



Gl 1 (E;k) 
Gf (E;k) 
Gf (E;k) 
Gf (E;k) 



E + E k cos 9 



E 2 

E- E k cos9 



El 



Ek 



= E k 



E 2 

sin 9 cos ( 



El 



i sm v sin < 



E 2 



'El 

i sin 9 sin < 



E 2 -El 



In order to include the effects of a boundary 
calized in real space, the bulk Green's function 
in real space via a Fourier transformation. 



G Q (E;y) 



B° fit ... 
B B" St 
: B B a fit 
B B° 

B J < 



B 



(A6a) 
(A6b) 
(A6c) 
(A6d) 

that is lo- 
is written 

(A7) 



for 



B l 



dk 
2^ 



G (E;k), B 



dk 
2^ 



e lk G {E;k). (A8) 



There are of course more non-zero matrix elements, cor- 
responding to the mixing of matrix elements separated 
by more than one lattice constant. However, for systems 
with nearest-neighbor interaction, the matrix elements 
contained in B° and B are the only ones needed to prove 
the existence of zero energy states. 



2. Green's function of the open system 

Next we write an expression V that deletes the cou- 
pling terms to create a system with edges. In a system 
with nearest-layer interactions, the only non-zero matrix 
elements are those between a single pair of neighboring 
layers: 



V = 







v b 



(A9) 



10 



Vb = — b ■ er. and for the Hamiltonian described in 
Eq. (|A4|I takes the form 



the two edges, the first two columns can be written as 



V b = 





-b r + b; 



-b r - h 




(A10) 



We are now ready to compute the Green's function for 
the open system. The Dyson equation gives an exact 
expression for the open Green's function G in terms of 
the periodic Green's function Go and the cuts V needed 
to take the system from one geometry to the other: 

G(E; y) = (I- G (E; y)V)- 1 G (E; y). (All) 



Substituting Eq. (|A"9|) and ([17]) into Eq. (fAlT]) gives the 
following condition for an edge state wavevector u: 



(1 - G V)u 



I - BVJ 



-B°Vf 



B°V b 



I - B^V b 



u = 0. (A12) 



The ellipses indicate the only non-zero sub-blocks: the 
first and last columns, and copies of the identity along 
the diagonal. 



3. Existence and spectrum of edge modes 



u = 



Lo = 



1 + (b r + bi)(c x - ic y ) 
(b r + &j)(c - c z ) 
(b r + bi)(a x - ia y ) 
(b r + bi)(a - a z ) . 

(b r - bi)(c + c z ) 
1 + (b r - bi)(c x + ic y ) 
(K - bi)(a + a z ) 
(b r - bi)(a x + ia y ) . 



(A14) 



(A15) 



Now we examine the conditions required for the system 
to have a left edge state, i.e. when the columns L\ and L2 
are linearly dependent. This is done by writing L\ and L2 
as a 4 x 2 matrix M and requiring that the determinant 
of any 2x2 sub-block vanish. This constraint takes on 
the following form for the sub-blocks of M created by the 
bottom two rows, the top two rows, and the second and 
fourth rows, respectively: 



a 2 , 



at), 



:(&* -&?)(<£+■_ 

1 + 2(b r c x - ibiCy) 

+ (& 2 -& 2 )( C 2 +c 2 -c 2 + C 2 ), 



(A16a) 
(A16b) 

= (b 2 r - b 2 )(c - c z )(a x + ia y ) 

- (b 2 r - b 2 )(c x + ic y )(a Q - a z ) - (b r + b l )(a - a z ). 

(A16c) 

The conditions are equivalent to B VfjuL = and (/ — 
BV b ^)u L = 0, which together require Det [B°] = 0. To 
see this, note that although the first condition is satisfied 
when either Det = or V b ul = 0, the second case 
cannot satisfy (/ — BV b )uL = 0. 



Since an edge state wavevector u has non-zero com- 
ponents in the sub-block corresponding to one of the 
edges, Eq. (| A12|) can be satisfied in two ways. There 
is an edge mode on the left edge when the first two 
columns of Eq. (|A12[) are linearly dependent, and on the 
right edge when the last two columns are linearly de- 
pendent. Recalling the expressions for the bulk Green's 
functions, Eq. (|A6|) . we introduce some notation for the 
open Green's function: 



a Q 



dk E 
27 £ 2 - E 2 k ' 
dk Ei. sin 6 cos ( 



2tt E 2 
dk Ek sin 9 sin < 



El 



2tt E 2 
dk E]~ cos 6 



El 



2tt E 2 



El 



(A13a) 
(A13b) 
(Al3c) 
(A13d) 



„ _ r dk Ak 
C - J 27 e 



and similar for c x ,c y , and c z . Let us 



E 2 -El> """"" iwi ^xi^yi 

first examine the conditions required for a left edge mode. 
Collapsing Eq. (|A12[) to the subspace corresponding to 



Constraints on hn and E 



We begin by showing that 
Det [B°] =a 2 - 



(A17) 



is zero if and only if hy encloses the origin and the 
edge state energy is given by E = ±|bj_|. For a ge- 
ometric view of the variables ao,a x , etc., note that 
Ek (sin 9 cos </>, sin 9 sin <f>) are the (x, y) coordinates of h|| , 
and E k cos 9 — b° ± -z. El is given by Ib'JJ 2 + |h|| | 2 , and 
E 2 — El < for a mid-gap state. We examine two cases, 
when hy encloses the origin and when it does not. 

Suppose the ellipse does not enclose the origin. Let 
h™ 1 " be the point on the ellipse hy closest to the origin. 
Because the ellipse is convex and does not contain the ori- 
gin, it must lie in the half of the plane (spanned by b r , b 4 ) 
for which r • h™ ln is positive. In other words, hy (k) ■ h™ 111 

is positive definite. Moreover, hy(/c)-h™ m > |h™ ln | 2 from 
its definition. The Cauchy-Schwarz inequality says that 



\(a x x + a y y) ■ hT n \ < \a x x + a y y\ |h]j 



(A18) 
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therefore 



to get: 



\a x x + a y y\ > 



> 



h " lln | 



dk hii • h : 



2-kEI 
dk \K 



E 2 



2ir El- E 2 ' 



(A19) 



a y y 



dk 



27T(|h|| 



(A22d) 



By the triangle inequality: |/dA;u(fc)| < fdk \u(k)\, we 
can put an upper bound 



This sets an upper bound to Det [B°] : 



ii-a 2 y <(E 2 -\b° ± \ 2 )l?-\h^\ 2 X 2 , 



(A20) 



where X = J 



dk 



Since |b°_| 



| h mi„|2 > E 2^ 



the 



2tt El~E 2 

expression is always negative and never zero. Hence no 
edge states can exist when the ellipse hy fails to enclose 
the origin. 

Now we consider the case when the ellipse hn encloses 
the origin and show that an edge state exists only when 
E 2 = \b\\ 2 . First, when E 2 = \b\\ 2 , the integrals a 
and a y are both zero. Because E 
and the determinant is zero. 

To see why a x = a y = in this case, note that the 
denominator of X becomes |h||| 2 , and we can express the 
conditions geometrically. Here we use Gauss's law in 2D 
to show that J dk is zero if and only if h|| encloses 
the origin. We can visualize the expression from an elec- 
trostatics point of view: for a charged ellipse in 2D with a 
charge distribution JdkS 2 (r— h||), the electric field (oc - 

in 2D) at the origin is given by J dk u 'U ■ Analogous to 

a uniformly charged circle, the charges are distributed 
such that the electric field is vanishing in the interior but 
nonzero in the exterior: 



a x x + a y y 



dk-, 



0. 



\a x x + a y y\ < e 



dk 



< e 
= eX 



27r(|h,|P- e a)|h,||2 
dk e 
M|h||| 2 -e 2 )|h||| 
dk 1 
M|h||| 2 - £ 2 ) 



(A23) 



We have used the fact that |h|| | > e to go from the second 
to third line. This implies a 2 
bound to Eq. (|A17|k a\ - a 2 



a 2 < e 2 X 2 and sets a lower 
2 -a 2 > Ofor \E\ > |b° I 



- a 



v 



Hence we have shown that Det [B°~\ 



if and only if 
E 2 = |bj_| and the ellipse hy encloses the origin. 

Turning now to the second constraint, Eq. (|A16b|) , we 
note that E = ±|bj_| implies that Cg — c 2 . In the follow- 
ing, we use the fact that c x = —ic y . Substituting this into 
Eq. (|A16b[) gives the constraint 1 + 2(b r c x — ibiC y ) = 0, 
which may be expressed as 



, dk ,, hn ■ b* 



= o. 



(A24) 



In order to satisfy Eq. (|A16a[) and Eq. (|Al6b[) simultane- 
ously, Eq. (|A24|) must be true if and only if hii encloses 
the origin. This can be shown for an arbitrary ellipse 
centered at bjj and with semimajor and semiminor axes 



(A21) given by 2b r and 2b* 



Now we show that Det [B°] = is not satisfied for 
any other value of E 2 . If E 2 < \b° ± \ 2 then \E\ < \b° ± \ = 
\E k sm6 sin <j>\. It follows that a 2 , - a 2 = {E 2 -\h\\ 2 )X 2 < 
0, hence there are no edge states. 



If E 2 
E 2 = |h 



> |b^| 2 , we let E 2 



|2_ 
E 2 > 0. 



such that E\ 



oq and a z can be written as: 



-a = 



dk E 
2^|h||| 2 -e 2 ' 
dk \b a ± \ 
2^|h||| 2 -e 2 ' 

which combine to give: 

al-al = (\E*-Kf)X*=eV. 

For a x and a y , we use the manipulation 

1 1 e 2 



(A22a) 
(A22b) 

(A22c) 



|h||| 2 -6 2 



l h lll 2 



e 2 )|h,| 



5. Sign of the energy 



The final constraint, Eq. (|A16c[) . determines the sign 
of the edge state energy. Note that the first two terms 
vanish because a x — a y = and c x = —ic y . Getting rid 
of the remaining term requires choosing a particular sign 
for the energy. The condition is satisfied with the choice 
a Q = a z . 

The constraints for the left edge can be summarized as 



la. E = +b°, ■ z. 



lb. 1 



0. 



2(b r c x - ibiCy) 

A similar calculation for the right edge, using the last 
two columns of Eq. (|A12j) . gives the following conditions: 

2a. E = -b°, ■ z, 



2b. 1 + 2(b r c x + ibiCy) = 0. 



Condition 2b. can be expressed as 1 — 2 J |p 



dk g — 'ik 



Nl 2 



0. Hence it is also satisfied when encloses the ori- 
gin, but with opposite orientation as for the left edge. 
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The conditions for edge modes are now equivalent to two 
statements: 

1. The function h|| must enclose the origin, 

2. The energy is given by E = ±b*| • z, with the sign 
determined by the orientation of loop hii . 

We have succeeded in analytically deriving the condi- 
tion for the 2L x 2L matrix representing the Hamiltonian 
for an open system to have zero eigenvalues. The power 
of the method lies in the fact that V has non-zero matrix 
elements only in the 4x4 subspace of electron operators 
at the two edges. Thus the effect of the boundary can 
be seen by examining the 4x4 subspace, which can be 
handled analytically. 



Hence we need to find an energy E and a set of A's all 
within the unit circle, such that their corresponding null 
vector u defined by (|B3|) sum to zero [Eq. (|Blb|) ]. 

The outline of the proof is as follows. First, we de- 
rive the particle-hole relationship between left edge and 
right edge modes. Second, we establish the algebraic re- 
lations between the A's and E of an edge state [Eq. (|B8[) ] . 
Third, we define complex functions L, L which represents 
the ellipse traced out by (3{e lk ). Next, we show (assum- 
ing an edge state exists) that the energy of an edge state 
is given by the displacement of the ellipse |b^|. We then 
proceed to prove Thm. la, the condition which governs 
the existence of an edge state (i.e. when all |A M | < 1). Fi- 
nally, we compute the edge states projectors [Eq. (IB15I) ] 
and determine the sign of the edge state energies, which 
completes the proof for Thm. lb. 



Appendix B: Proof by transfer matrices 

In this section, we prove Theorem 1 for a hard-edge 
at a fixed k|| . We begin by defining a new function (3 
(which is like a complex extension of h) , and the form of 
our edge states. 

In the Hamiltonian (fTJ), b • T, b° ■ T, b* ■ T arc the hop- 
ping matrices from the previous, same, and next layers 
respectively. We consider a semi-infinite system where 
the layers are labelled by positive integers with n = 1 
denoting the layer at the surface. Hence we ignore the 
terms V&i^q an< ^ ^o^i m considering our semi-infinite 
system. Consequently, an excitation ip' = ^ n ^^ipn of 
the Hamiltonian (fTJ) satisfies the following properties: 

b • Ityn-i + b° • Tip n + b* • TVwi = Eip n n > 1, 

(Bla) 

V>o = 0. (Bib) 

An edge state ip is one for which i\) n is exponentially 
decreasing as a function of n. Due to the translational 
invariance in the bulk, we use the ansatz 

^n = $>M*2 (B2) 

in our solutions4^i An edge solution requires that < 
1 for all fi and that ^ M = 0. 
For each decaying mode (A, u) we have^ 

A [/3(A) ■ T - E] u = 0, (B3) 

where we have defined 

/3(A) = A -1 b + b°+Ab*. (B4) 

Notice that when A has unit modulus, we recover the 
Bloch equation for the bulk (propagating) modes of the 
system: 



1. Relation between left-right boundaries 

Lemma. For every left edge state with energy E, there 
is a corresponding right edge state with energy —E, and 
vice versa. 

The recursion relation (|Bla[) and boundary condition 
(|Blb[l describe a semi-infinite system with a "left" edge. 
We can write a similar system for the "right" edge sim- 
ply changing the condition in (|Bla[) to n < —1 xor swap- 
ping b with b* in the equation. The first transformation 
amounts to finding a set of A's outside the unit circle, 
i.e. \Xfil > 1, such that the edge wavefunction ip n de- 
cays with decreasing n. The second transformation is 
equivalent to taking the complex conjugate of b, and the 
two transformation are the same owing to the fact that 
/3*(A)=/3(1/A*). 

When the gamma matrices are the Pauli matrices, 
there is a charge conjugation operator C = —ia y K, which 
takes b cr > — b* er, where K. is the complex conjugation 
operator. Since C/3(A) • erC -1 = -/3*(A*) • <x, C turns a 
left edge state with energy E into a right edge state with 
energy —E and vice versa. 

When the Dirac matrices are larger then 2x2, such a 
C operator still exists, as it is always possible to find an 
antiunitary operator which flips the sign of three of the 
gamma matrices. This operator C will depend on what 
b and b° are, meaning that C is a function of ky , making 
it a non-local operator. However, the conclusion remains 
the same. 

For the remainder of the proof, we will only be focused 
on left edge states. 

2. Algebraic relation between A a ,Af, and E 

For any E, there are (at most) four possible A's sat- 
isfying (|B3[) . evident from squaring A/3(A) ■ T to get the 
quartic equation 



13 (e^) = h(k ± ). 



(B5) 



A/3(A) • A/3(A) - E 2 \ 2 = 0. (B6) 
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Note that if A is a root to this equation, then so is 1/A*. 
Hence there can be at most two solutions of for A within 
the unit circle, which we call A a and A;,. 

The edge wavefunction takes the form ip n = u A™ + 
w&A™ with u a = — u b to satisfy the hard-edge boundary 
condition (|Blb[) . where the coefficient u a is a (right) null 
vector of the matrix X a ((3(X a )-T— E) and similarly for u b . 
It follows that the matrices X a ((3(X a )~- E) and X b ((3(X b ) — 
E) must share a non-zero null vector, or equivalcntly, any 
linear combinations of the two matrices must be non- 
invertible. In other words, an edge state at energy E 
exists only if the following conditions^ are satisfied: 



Det [X a (fi(X a )-T-E)]=Q, 
Det [X b (fi(X b ) • T - E)} = 0, 
Det [c a X a ((3(X a ) ■T-E) + c b X b {fi{X b ) 



(B7a) 
(B7b) 

T-E)] = 0, 
(B7c) 



quadratic polynomials in A: 

L(X) = qX + w+pX^ 1 , 
L(X) = p*X + w*+ q*X _1 . 

w = bjj • (vi + iv2) and q = 



(BIO) 



where p = b • (vi + IV2) 
b* ■ (vi + iv 2 ). 

It is straightforward to show from (|B9j) that fiu (A ) • 
/3||(A&) = \(L a L b +L a L b ), where L(X a ) is abbreviated as 
L a , etc. Equations (|B8|) become 



X a L a L a , 



X 2 a (E 2 - A 2 ) = XI 
X 2 b (E 2 - A 2 ) = A 2 L 6 Z 



X a X b (E 2 - A 2 ) = h^(L a L b 



L b L a ), 



(Blla) 
(Bllb) 
(Bile) 



where A 



|b°J 



for arbitrary c a , c b and for | A a | , | X b | < 1 . The converse 
statement is also true, as (|B7a[) and (|B7b[) implies that 
the ranks of the matrices X a , b (fi(X a . b ) • T — E) are at 
most half their dimension (a property of gamma matri- 
ces). The last equation (|B7cj) means that the two ma- 
trices must share a right null vector or a left null vector. 
In the former case we have a left edge state at energy 
E, and in the latter case we have a right edge state at 
energy E, which by our lemma implies a left edge state 
at energy — E. 

We may rewrite the equations in a more useful form: 



A 2 £ 2 = A 2 /3(A a )-/3(A a ), 
A 2 £ 2 = A 2 /3(A b )-/3(A b ), 
X a X b E 2 = X a X b f3(X a ) ■ f3(X b ). 



(B8a) 
(B8b) 
(B8c) 



3. Introducing functions L, L 



bj_, we de- 



Similar to b° = bjj + bj_ and h = ii|| 
compose fi into components parallel and perpendicular 
to the 2D plane (ID if the ellipse is degenerate) spanned 
by {b r , b 1 }: /3(A) = fi\\ (A) + b°_ Keep in mind that while 
bj_ is a real vector, fin (A) = A _1 b + bjj + Ab* is generally 



a complex vector, unless A 



■4k 



has unit modulus. 



We want to find complex functions which trace out 
the same ellipse as f3\\(e lk ) in the complex plane. We 
first choose two real orthogonal unit vectors Vi , v 2 as a 
coordinate basisof the 2D plane. Let 



L(X) 
L(X) 



0|| (A) 
0|| (A) 



vi +ij3||(A) 
vi-i/3,|(A) 



V2, 
V 2 . 



(B9) 



The loci L(e lk ) and L(e lk ) both trace out the ellipse in 
the complex plane identical to @\\(e ik ), but with different 
orientations. In general, L(X) and i(A) are not conjugate 
pairs unless A lies on the unit circle. Expanding their 
definitions, we can see that both AL(A) and AL(A) are 



4. Edge state energy 

In this portion, we show that the existence of an edge 
state requires: E = ±A, |A a ,b| < 1, and either L a = L b = 
or L a = L b = 0. The converse statement is trivially 
true by inspecting Eq. (jBllj) . Hence if an edge state of 
the semi-infinite chain exists, we show here that it must 
have energy E = ±|bj_|. (First half of Thm. lb.) 

Assuming that neither A a or X b arc zero, then 



E 2 
E 2 



A 2 
A 2 



L a L a , 
L b L b , 



E 2 -A 2 = \{L a L b + L b L a ). 



which we can combine to get 



(L a - L b ){L a - L b ) = 0. 



(B12a) 
(B12b) 
(B12c) 

(B13) 



Equations (|B12|) and (|B13|) are simply reformulations 
of the recursion relation (|Bla[) and boundary condition 
(|Blb|) . We now proceed to show that E = ±A. 

Proof by contradiction Suppose that E 2 — A 2 7^ 0. 
Then L a ,L bl L a ,L b are all non-zero. Eq. (|B13|) implies 
that L a = L b or L a = L b . Eq. (|B12a|l and (|B12b^ to- 
gether means that one equality implies the other, and 
hence L a = L b and L a = L b are both true. Here we have 
two polynomials with roots A a and X b : 

qX 2 + (w- L a )X+p = 0, 
p*A 2 + (w* -L a )X + q* = 0, 

which means that p/q = q* /p* = X a X b . However, since 
|A Q Af,| < 1, we have a contradiction as \p/q\ cannot be less 
than one and greater than one at the same time. Hence 
we have shown, should an edge state exist, it must have 
energy E = ±A. Now (|B12a[) tells us that either L a or 
L a is zero which combined with (|B13[) leads to the desired 
result. 

If one of A a ,Ab is zero, say A& = (which happens 
when the ellipse is a circle), then the expressions simplify 
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as XbLb = p and XbLb = q*. From (lBllb[) . either XbLb or 
AbLb is zero. From (|Bllc[) we have either L a = Lb = 
or L a = Lb = 0, and in either case, (|Blla|) implies that 
E 2 - A 2 = 0. 

If both A a and A& are zero (which happens when the 
circle is centered on the origin), then A = must be a 
double root to the polynomial (|B6I) . In this case, we have 
a flat band, and it is much easier to refer back to (|B1[) 
and solve the system directly. One easily finds that the 
statement about edge states holds. 



5. Existence of edge states 

Finally, we use the fact that |A |, |A&| < 1 to determine 
when an edge mode is present. Recall that there are four 
zeroes (and two poles) to the equation L(X)L(X) = 0, 
and that at most two of the roots have modulus less than 
one. An edge mode exists if either L(X) or L(X) has both 
roots A within the unit circle. (No edge mode exists if 
each function L,L has one root within the unit circle.) 

We can compute the number of zeroes (inside the 
unit circle) of the function L by the contour integral 

27tI $ L L(z) ^ Z a l° n § the unit circle, which computes the 
number of zeroes minus number of poles within the unit 
circle. As L{\) has one pole (at A = 0), L(X) has two 
zeroes if and only if the ellipse L{e lk ) wraps around the 
origin counterclockwise, leading to an edge state. On 
the other hand, if L wraps around the origin clockwise, 
then L(e ik ) wraps around the origin counterclockwise, 
and there are two zeroes for L(X) within the unit circle 
which also leads to an edge state. In the case where the 
ellipse /3|| does not wrap the origin, then neither L or 
L has two roots within the unit circle, and an exponen- 
tially decaying solution to the semi-infinite system does 
not exist. 

This completes the proof for Theorem la, which relates 
the presence of edge states to the properties of the ellipse 
/3(e ife ) = h(fc). 



6. Sign of edge state energy 



b°_| > 0. This 



In this section, we determine whether the left edge 
state energy is +A or —A, where A 
will complete the final statement of Thm. lb 

Define vj_ as the unit vector parallel to (assume 
A > 0). Recall that Vi,V2 are unit vectors used in 
the definitions of L, L, and so the three unit vectors are 
mutually orthogonal. Define the corresponding gamma 
matrices T x = vi • r,r* = v 2 ■ I\ Y 1 - = v_l ■ I\_which 
pairwise anticommute. As (3\\ ■ Vi = ^(L + L) and 
/3 |r v 2 = i(L-Z), 



A(/3(A) ■T-E) = X(AY ± - E) 



\L{\) 



(Y x - iY y ) 



AL(A) 



(Y x + iY y ). (B14) 



The first term annihilates the projector ^(1 + ^Y- 1 ) while 
the second and third terms annihilate ^(1 + iY x Y y ) and 
i(l + iY x Yy) respectively. 

By inspection, the projector P = J2 uut 011 to the edge 
state depends on whether L or L has two zeroes inside 
the unit circle, as well as the energy E: 



P± = 



} i(l±r ± )(l-iFT*') L a = L b = Q 
\(l±Y ± )(l + iY x Yy) L a = L b = 



(B15) 



where P± projects on to the edge states with energy E = 
±A. 

When the irreducible representations of Y arc 2x2, the 
product iY x Y v must equal either Y 1 - or — T^, hence cither 
P + or P_ must be zero. This implies that an edge state 
occurs only at A or —A but not both. For example, when 
T x T y = iT ± and £ = o, then P+ = |(l+r- L ),P_ = and 
there is exactly one edge state at energy A. In general, 
we can determine E via the orientation of the ellipse h: 



E = b°- 



h r x h l 
\h r x 



(B16) 



where the cross product is defined from the commutation 
algebra b r x b' = Tr{cr [b r • ctM ■ a}}. This result 
gives us Theorem lb. 

When the irreducible representations of Y are 2 m x 2™ 
with m > 2, there are edge states at both A and — A. As 
P+ + P_ = i(l =F iY x Yy), there are a total of 2™- 1 left 
edge states. Notice that the projectors P + and P_ are 
related by the similarity transformation P + = Y W P_Y W , 
where Y w anticommutes with Y x,v,± , and so there must 
be 2 m_2 edge states at each energy. 



7. Effective surface Hamiltonian 

We can use the edge state projector to construct the 
effective surface Hamiltonian. The operator 

P s = P+ + P- = |(1 - iY x Y y ) (B17) 

projects on to the edge states (assuming vi and V2 are 
chosen such that L = 0). Notice that it only depends on 
b ■ r, the nearest-neighbor coupling and not the on-site 
potential: 
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i / fb r • r b* • n 

1 - tY x YV) = - 1 - i l ,. ' . J 



in 



2|b r x b 

[b-r,b* -r] 



2 \ — (Area of ellipse) 
The effective surface Hamiltonian is 



(B18) 



H s = EP+ - EP- + £co(l - P s ) 

= P s (b° L -Y) + E oc (l-P s ), (B19) 

where E^ — > oo such that the low energy theory de- 
scribes the surface states. 
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Appendix C: Proof for continuum Hamiltonians 



The proof of Theorem 3 is very similar to the proof 
of Theorem 1 using transfer matrices, and so we present 
here a condensed version of the proof. 

Consider a semi- infinite system with x > 0, and cither 
■0(0) = (Dirichlet) or i/>'(0) = (Neumann) boundary 
condition. The momentum parallel to the surface is a 
good quantum number, and so we fix p|| to get an ef- 
fective ID problem. We seek a solution of the form^S 
ip(x) = ^^c^^'u^, with Imft M > 0. Each pair (k,u) 
satisfies: 

(C*° + C x k + C 2 k 2 - E)u = (h(re) -T-E)u = Q. (CI) 

Squaring h • T gives us the quartic equation 

h(«) • h(«) - E 2 = (C2) 

with real coefficients. Hence if k is a root, then k* is 
also a root. Once again, we have at most two solutions 
for k in the upper half of the complex plane, and so 
the wavefunction must take the form %p(x) = u a e lKaX + 
u b e lKbX . With either Dirichlet (u a = —Ub) or Neumann 
{n a u a = —KtUb) boundary condition, we have u a oc Ub 
and so h(re a ) - T — E and h(nb) T — E share a null vector. 
By a similar argument to that in Sec. IB 21 the existence 
of an edge state is equivalent to 

Det [c a (h(K a ) -T-E) + c b (h(K b ) ■ T — E)] = 0, (C3) 

for all c a ,Cb and ImK a ,ImK(, > 0. Equivalcntly, we have 

E 2 = h(« a ) • h(« a ) = h(Kb) ■ h(/t 6 ) = h(n a ) ■ h(n b ). 

(C4) 

We introduce the functions L(k),L(k): 



L(k) = h|| («) • vi + ih|| (re) ■ v 2 , 
L(k) = hii (re) • vi — ihn (re) ■ V2, 



(C5) 



where vi , v 2 form an orthonormal coordinate basis in the 
plane spanned by c 1 and c 2 . L(re),Z(re) are quadratic 
polynomials in k: 

L(k) = (c° + ic° y ) + (4 + »cJ)k + (c 2 + lC 2 ) K 2 , 
L(«) = (c° - t c°) + (4 - ici)« + (c 2 - lC 2 ) K 2 . 



When k = p is real, £(p) and L(p) trace out the parabola 
h|| (p) in the complex plane with opposite orientations. 
Using the relation hii(« a ) ■ h|i («b) = k(L a Lb + L a L b ), 
where L(K a ) is abbreviated as L a , etc., Eq. (|C4|) become 



E 2 - A 2 = i a Z a , 
iJ 2 - A 2 = L b Z b , 
E — A = \{L a L b + LbL a ), 



where A = | c j_ | . The equations combine to get 
(L a -L b ){L a -L b ) = 0, 



(C7a) 
(C7b) 
(C7c) 

(C8) 



Equations (|C7ajl . (|C7bj) and (fC8|) together with Im n a ^ b > 
are true if and only if an edge state exists at energy ±E. 

We construct a proof by contradiction showing that 
E = ±A 2 . Suppose E 2 £ A 2 , the n all of L a ,L b ,L a ,L b 
are nonzero. Equating (|C7a[) and (|C7b[) gives L a /L b = 
La/Lb, and combining with (|C8[) implies L a = Lb and 
L a = Lb- Hence the polynomials L(k) — L a and L(k) — X a 
have identical roots («;„ and Kb). The sum of the roots 
K a + K b must lie in the upper half plane, and it is equal 
to (4 + ic\)l(c 2 x + ic 2 y ) and {c\ - ic y )l(c 2 x - ic 2 y ) from 



C6|) . This leads to a contradiction as the expressions 
are complex conjugate pairs. Therefore, an edge state 
requires E = ±A and either L a = L b = or L a = Lb = 
0. 

Finally we impose the condition that ImK Q , lmn b > 0. 
There are no poles in the function L(k), and so the 
number of zeroes in the upper half plane is given by 

2^7 § L{z) d z i integrated along the real line from — R to 
R and closed on the upper half plane Re 10 for < 9 < ir, 
where R is taken to infinity. Assuming that c 2 is non- 
zero and so L is a quadratic function of z, the contour of 
L(z) for z — Re 10 always wraps the origin by 2ir radians 
counterclockwise. 

Hence L(k) has two roots in the upper half plane if and 
only if the parabola L(p) wounds around the origin coun- 
terclockwise. Similarly, L(k) has zero roots (so L(k) has 
two roots) in the upper half plane if the parabola wounds 
around the origin clockwise. An edge state exists in both 
these cases, which occur when the origin lies in the con- 
cave side of hj| (p). If the origin is not in the concave side 
of the parabola h|| (p), then L(k) and L(k) only have one 
root in the upper half plane and the system has no edge 
states. This completes the proof for Theorem 3a. 

To determine the sign of the edge states, we construct 
the projectors for E = ±A. The projectors in the con- 
tinuum case is identical to that of the lattice case (|B15|) , 
hence by the same argument used for Thm. lb, we can 
prove Theorem 3b. 



E 



C 1 X c z 



(C9) 



In addition, one may also derive the effective surface 
Hamiltonian: 

H s = EP+ - EP- + -Eoo(l - Ps) 

= P s (c° L -T) + E oc (l-P s ), (C10) 
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A solution of the form u^nA'^ 1 is also possible, when A M 
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